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The properties of hybrid stars formed by hadronic and quark matter in /3-equilibrium are described 
by appropriate equations of state (EoS) in the framework of the quark meson coupling (QMC) model. 
In the present work we include the possibility of trapped neutrinos in the equation of state and obtain 
the properties of the related hybrid stars. We use the quark meson coupling model for the hadron 
matter and two possibilities for the quark matter phase, namely, the unpaired quark phase and the 
color-flavor locked phase. The differences are discussed and a comparison with other relativistic EoS 
■ is done. 
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During the early stage of a proto-neutron star neutrinos get trapped in it when their mean-free path is smaller than 
I the star radius. The presence of neutrinos generally gives rise to a stiffer EoS. This may have important consequences 
on the evolution of the star, namely, it can happen that during the cooling process the star decays into a low-mass 
black hole [lj . It is worth mentioning that over 20 years ago the importance of neutrino trapping was already pointed 
out in , where the authors claimed that another important effect of including trapped neutrinos is that the collapse 
' 1 is gentler in its presence than it would be without it. 
, In the present paper we are interested in buildingthe neutrino trapped EoS for mixed matter of quark and hadron 
phases. We employ the QMC model (QMC) 0. [113 in order to describe the hadron phase. For the quarkphase we 
| use two distinct models, the unpaired quark model (UQM), which is given by the simple MIT bag model |(j and the 
t-H , color-flavor locked phase (CFL) |7], |8( in which quarks are paired near the Fermi surface forming a superconducting 
phase [jj. 

In a previous work [To| we have used the same formalism in order to study the properties of hybrid stars at T = 
MeV. In this work, we verify the effect of including trapped neutrinos in the same spirit as done in where we 
have seen that the EoS changes considerably and the mixed phase appears at higher energy densities than in the EoS 
— . ■ built without the inclusion of neutrinos in accordance with what was seen in other works, as in [lj for instance. We 
have also verified that trapping keeps the electron population high so that dense matter contains more protons (and 
depending on the parametrization used, other positively charged particles) than matter without neutrinos. This fact 
O ' was also discussed in |12|. 

Proto-neutron stars with a certain baryonic mass at birth keep this mass during its evolution until the final neutrino 
free star because most of its matter is accreted in the early stages after birth. Hence, stars with trapped neutrinos 
and a baryonic mass larger than the corresponding ones after deleptonization collapse to a black hole. In 11] we have 
described the hybrid star within a non-linear Walecka model (NLWM) ^] including the baryonic octet and a phase 
transition to a quark phase. Within this description the maximum baryonic masses supported by neutrino trapped 
EoSs are larger than the corresponding ones found in neutrino free EoSs, a fact which occurs in other EoS which 
include the strangeness degree of freedom P, • I n what follows we investigate whether this behavior is also present 
in the framework of the QMC model. Because of its importance in understanding the evolution process in a star, we 
calculate the baryonic masses for the hybrid stars studied in the present work. We also compare the properties of the 
stars obtained within the UQM model and the CFL model. 

We perform all the calculations at T = MeV although the neutrino trapped phase occurs for temperatures in 
the interior of the star which can vary between 20-40 MeV 0, ^J. However, it was shown in llj that the effect of 
trapping is much stronger than the finite temperature effect. Therefore, we believe that the main conclusions drawn 
in the present work are still valid for a finite temperature calculation. On the other hand, the calculation with the 
CFL phase will give us only an upper limit, since at finite temperature there is a phase transition from the color 
superconducting state to a normal phase, which, according to |l6j |. is a second order transition with a BCS critical 
temperature T c ~ 0.57A, A being the gap parameter. 

Since most of the analytical calculations and formulae used in the present work have already been given in [Toj 
they will be omitted here. As mentioned above, we have used the QMC model with the inclusion of hyperons for the 
hadron phase. In this model, the nucleon in nuclear medium is assumed to be a static spherical MIT bag in which 
quarks interact with the scalar and vector fields, a, lo and p and these fields are treated as classical fields in the mean 
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field approximation. The quark field, ip q (x), inside the bag then satisfies the equation of motion: 

i p - (m° - gl a) - fl« wt° + \g*T tPo3 ^{x) = 

q = u,d,s, (f) 

where m° is the current quark mass, and g%, g% and denote the quark-meson coupling constants. After enforcing 
the boundary condition at the bag surface, the transcendental equation for the ground state solution of the quark (in 



s-state) is jo(x q ) = f3 q ji(x q ) which determines the bag eigenfrequency x q , where (3 q = ■J (fl q — Rbttl*) / (Q q + Rsm*), 
with Q q = (x q + R 2 m* q Y^ 2 'i m q — m q ~ 9a a o, is the effective quark mass. The energy of the nucleon bag is 

M B = 3^r- — jf- + ^ttR b Bb, where Bb is the bag constant and Zb parameterizes the sum of the center-of-mass 
motion and the gluonic corrections. The bag radius, Rb, is then obtained through the stability condition for the 
bag. An interesting fact related to the QMC model is that the bag volume changes in the medium through the mean 
value of the c-field. This also implies that the bag eigenvalues are also modified. The onset of hyperons depends on 
the conditions of chemical equilibrium and charge neutrality discussed below and also on the meson-hyperon coupling 
constants for which we have chosen the hyperon coupling constants constrained by the binding of the A hyperon in 
nuclear matter, hyper-nuclear levels and neutron star masses (x j r = 0.7 and x u = x p — 0.783) and assumed that the 
couplings to the £ and S are equal to those of the A hyperon [13j . The leptons either in the hadron or in the quark 
phase are considered as free Fermi gases and they do not couple with the hadrons, with the mesons or with the quarks. 

The condition of chemical equilibrium is imposed through the two independent chemical potentials for neutrons fi n 
and electrons fi e and it implies that the chemical potential of baryon Bi is /is 4 = Qf ' fi n — Qlli e , where Q\ and Qf are, 
respectively, the electric and baryonic charge of baryon or quark i. Charge neutrality implies J2 B QiPBi + J2i UPl = 
where qi stands for the electric charges of leptons. If neutrino trapping is imposed to the system, the beta equilibrium 
condition is altered to fiB t = Qf fi n — Q1(fi e — fi Ve ). In this work we have not included trapped muon neutrinos. 
Because of the imposition of trapping the total leptonic number is conserved, i.e., Yj, = Y e + Y Ve = 0.4. 

For the quark phase we consider two models. First of all, we take the quark matter EoS as in 6] in which u,d and 
s quark degrees of freedom are included in addition to electrons. Up and down quark masses are set to zero and the 
strange quark mass is taken to be either 150 or 200 MeV so that we are able check the effect of the s-quark mass. In 
chemical equilibrium fid = fi s = fi u + fi e . In terms of neutron and electric charge chemical potentials fi n and fi e , one 
has fjL u = |/Lt n — |/i e , fid = \pn + \^e: Ps = \pn + \fie- In the energy density for the quark matter EoS a term 
+B and in the pressure a term — B are inserted. This term is responsible for the simulation of confinement. For the 
Bag model, we have taken B 1 / 4 =190 and 200 MeV. 

In the EoS taking into consideration a CFL quark paired phase, the quark matter is treated as a Fermi sea of 
free quarks with an additional contribution to the pressure arising from the formation of the CFL condensates. The 
density of the three types of quarks are identical and the electron density is zero, as shown in ML The expressions for 
the energy density and pressure depend on a gap parameter A which is taken to be 100 MeV 0. 

Once the hadron and quark phases are well established, we have to construct the mixed phase, imposing charge 
neutrality globally, XP? P + (1~ X)P^ P + P l c = 0, where p l c p is the charge density of the phase i, x is the volume fraction 
occupied by the quark phase and p l c is the electric charge density of leptons. We consider a uniform background of 
leptons in the mixed phase since Coulomb interaction has not been taken into account. According to the Gibbs 
conditions for phase coexistence, the baryon chemical potentials, temperatures and pressures have to be identical in 
both phases, i.e., fi H p, n = PQP.n = Pn, Php, & = Pqp,e = Pe, T HP = T QP , P HP {fi ni fi e T) = P QP {fi ni fi e ,T), 
reflecting the needs of chemical, thermal and mechanical equilibrium, respectively. 

In fig. ^ the EoSs obtained with both quark models are displayed for different B values and two strange quark 
masses with neutrino trapping (Yl = 0.4). For the sake of comparison we have also plotted one EoS with no neutrinos 
(Yu e — 0), and included the EoS obtained within a NLWM formalism for the hadronic phase [ll| and an UQM with 
m s = 150 MeV and f? 1 / 4 = 190MeV, with and without neutrino trapping. As already discussed in 0, the EoSs 
are harder if neutrino trapping is imposed, independently of the model used. A larger s-quark mass and a larger B 
parameter make the quark EoSs harder in the mixed phase, a fact that manifests itself on the maximum mass stellar 
configuration. The main differences between the QMC formalism and the NLWM are: a) the NLWM EoS is harder 
at low densities and softer at intermediate densities due to the presence of hyperons; b) the transition to a pure quark 
phase occurs at lower densities in the NWLM. This behavior has consequences on the properties of the corresponding 
families of stars. 

In order to better understand the importance of the neutrinos when neutrino trapping is imposed, in fig. [21 the 
fraction of neutrinos is shown. The behavior encountered for the neutrino fraction if the UQM is used resembles the 
one shown in 11]: the population of neutrinos decreases in the hadron phase and, contrary to |ll|. only increases 
in the mixed and quark phases. The highest yields are of the order of 0.16. In this model, for smaller values of the 
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TABLE I: Star properties 
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strange quark mass, the neutrino population at high densities is greater. This occurs because a smaller s-quark mass 
gives larger s-quark densities and therefore, a smaller electron fraction. A fixed lepton fraction then implies a larger 
v e fraction. 

If the CFL is chosen for the quark phase, the population of neutrino is higher in the mixed and quark phases. This 
is due to the fact that in the CFL phase no electrons are present since the number of u, d and s quarks are equal, 
therefore the lepton fraction is kept constant only by the presence of neutrinos. This implies a much greater neutrino 
flux during dcleptonization. However, in a self-consistent calculation at finite temperature we do not expect such a 
strong effect since pairing will be weaker. 

The amount of neutrinos depends on the fraction of charged hyperons and quarks present in each phase, which are 
determined by the models used and consequently by the bag pressure and the strange quark mass. In the present 
approach with the CFL phase or UQM for i? 1 / 4 = 190 MeV the onset of hyperons, if any, occurs for p > 10po- 
With the UQM and B 1 / 4 — 200 MeV the hyperon onset occurs at ~ 4/jn but the charged hyperon fraction, namely 
£~, is never larger than 0.004. Hybrid neutron star profiles can be obtained from all the EoS studied by solving 
the Tolman-Oppenheimer-Volkoff equations. Even at finite temperature the conditions of hydrostatic equilibrium are 
nearly fulfilled 15]. In Table I we show the values obtained for the maximum gravitational and baryonic masses and 
radii of neutron stars as function of the central density for the EoSs studied in this work. For a fixed bag constant, the 
stellar and baryonic masses of the most massive stable stars are higher for higher strange quark masses in both models. 
For the UQM, these critical masses are always higher than for the CFL model, because it also corresponds to harder 
EoS. The radii and central energy density depend on the model and on the strange quark mass. The radii values are 
larger for the UQM model and the central energy density are larger for the CFL model, again due to the fact that 
the UQM EoSs are harder. Comparing the results of table I with the ones presented in [10(, where neutrino trapping 
was not considered, one can see that the inclusion of trapping makes the gravitational masses reasonably higher. The 
same conclusion was drawn in |ll| where the calculations were performed with different relativistic models. In the 
Table I we have also included properties of maximum mass stars obtained within the NLWM for the hadronic phase 
[ll| . Two conclusions are in order: the maximum baryonic masses obtained within the NLWM are larger and the 
difference of maximum masses for trapped and untrapped matter is smaller for the QMC (~ 0.2Af Q ) than for the 
NLWM (~ 0.4Af Q ). This means that the number of stars that would decay into a blackhole is much smaller in the 
QMC model and is probably due to the fact that no hyperons are formed in the interior of stars obtained with QMC 
for m s — 150 MeV and i? 1 / 4 = 190 MeV contrary to the NLWM case. If the quark phase is a CFL state the baryonic 
mass difference between the neutrino rich stars and neutrino poor is greater than in the UQM, (~ 0.35M©). This is 
understood because the greater flux of neutrinos carries out more energy. However for a finite temperature calculation 
we expect a smaller effect. 

A similar analysis was done in |14|. where the authors used a derivative coupling model with hyperons for the 
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hadron phase and the Bag model for the quark matter. They did not obtain any mixed phase for bag values larger 
than B 1 / 4 = 190 MeV in contrast with the present work. Moreover, the maximum masses shown in [l4j (~ 1.6M Q ) 
and the differences between maximum masses in neutrino rich and neutrino poor stelar matter are lower than in our 
calculations. 

In fig. |3]we display the baryonic masses versus the gravitational masses for both models. It is seen that neutrino 
trapped EoSs give rise to greater gravitational masses for the same baryonic mass. The mass difference reflects the 
binding energy released during the deleptonization and cooling stage . The maximum baryonic mass of the neutrino 
rich EoS are larger than the neutrino poor. This will lead to a blackhole formation during the leptonization period 
for stars with baryonic masses greater than the maximum baryonic mass of the neutrino poor EoS. This behavior has 
been encountered in other EoS which include strange matter, namely hyperon and/or quark matter |lllll2j . 

In summary, we have investigated the effects of neutrino trapping in the properties of neutron stars within the QMC 
framework, including the possibility of hyperon formation and a transition to an unpaired quark-phase or a CFL phase. 
We have concluded that within the QMC model with hyperons for the hadronic matter, either the hyperons only occur 
at very high densities or in very small amounts at lower densities (e.g. UQM). Another important point is that the 
maximum mass of a neutrino rich neutron star decreases after neutrino diffusion leading to the formation of a low 
mass black-hole. This mass reduction is smaller for a quark phase described within an unpaired quark phase than a 
CFL phase. If the quark phase is in a CFL state a large fraction of neutrinos is expected in the mixed and quark 
phases which will carry away more energy as they diffuse out. We also point out that the amount of neutrinos present 
in the CFL phase is almost the double in comparison with the amount found in the UQM phase. At finite temperature 
the effect will not be so strong and a self-consistent finite temperature calculation should be performed. We have 
also seen that the mass reduction of the maximum mass stars during to neutrino diffusion is smaller within a QMC 
formalism than in a NLWM formalism. 

This work was partially supported by CNPq (Brazil), CAPES (Brazil), GRICES (Portugal) under project 100/03 
and FEDER/FCT (Portugal) under the project POCTI/35308/FIS/2000. 
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FIG. 1: Equation of state obtained with the QMC model plus (a) UQM (b) CFL. 
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FIG. 2: Neutrino fraction for the EoS obtained with the QMC model plus (a) UQM (b) CFL. 
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FIG. 3: Baryonic mass versus gravitational mass obtained from the EoS with QMC plus (a) UQM (b) CFL. 



